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Observation of kink waves in solar spicules 
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Abstract. Height series of Ka spectra in solar limb spicules obtained with the 53 cm coronagraph of the Abastumani 
Astrophysical Observatory are analyzed. Each height series covered 8 different heights beginning at 3800 km above the pho- 
tosphere. The spatial difference between neighboring heights was 1", consequently ~ 3800 - 8700 km distance above the 
photosphere has been covered. The total time duration of each height series was 7 s. We found that nearly 20% of measured 
height series show a periodic spatial distribution of Doppler velocities. We suggest that this spatial periodicity in Doppler ve- 
locity is caused by propagating kink waves in spicules. The wave length is found to be ~ 3500 km. However the wave length 
tends to be ~ 1000 km at the photosphere due to the height variation of the kink speed. This probably indicates to a granular 
origin for the waves. The period of waves is estimated to be in the range of 35-70 s. These waves may carry photospheric energy 
into the corona, therefore can be of importance in coronal heating. 
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1. Introduction 

The heating of the upper chromosphere and corona is still 
an unsolved problem in solar physics. It is clear that the 
energy source supporting the high temperature lays in the 
highly dynamical and dense photosphere. There the mechan- 
ical energy of photospheric motions can be guided upwards 
by structured magnetic fields in the form of waves or elec- 
tric currents. Therefore the energy transport by magnetohy- 
drodynamic (MHD) waves throughout the solar atmosphere 
is one of the key process towards the solution of the heat- 
ing problem (Roberts 2004 1. Observations of oscillatory phe- 
nomena in the solar atmosphere have increased dramatically 
in the last few years through observations from SOHO (Solar 
and Heliospheric Observatory) and TRACE (Transition Region 
and Coronal Explorer). The space missions uncovered the rich 
spectrum of MHD oscillations in the transition region and 
corona (Doyle et al. J-999 Nakariakov et al. J^999 Ofman et 
al.[2000 Banerjee et al. 120011 O'Shea et al. 120011 De Moortel 
et al. i200 4). The coronal waves are either generated in situ or 
they penetrate from the photosphere. Hence the observation of 
oscillatory phenomena in chromospheric spectral lines is of vi- 
tal importance. 

Most of the chromospheric radiation in quiet Sun regions 
comes from spicules, which are jet-like chromospheric struc- 
tures observed at the solar limb mainly in Ha line. Spicules 
often show the group behaviour and probably are concentrated 



between supergranule cells (see e.g. reviews of Beckers [19721 
and Sterling 2000 1. Therefore spicules probably are formed 
in regions of magnetic field concentration and consequently 
MHD wave propagation in the solar atmosphere may be traced 
through their dynamics. Oscillations in spicules with ~5 minute 
period have been detected by ground based (Kulidzanishvili 
& Zhugzhda 1 983 ) and recently by space observations (De 



Pontieu et al. 120031 Xia et al. 2005 1. On the-other-hand, oscil- 
lations in spicules with shorter period have been reported more 
than 30 years ago by Nikolsky & Platova ( 1971 ). They found 
that spicules oscillate along the limb with a characteristic pe- 
riod of about 1 min. The oscillation, which they reported, was 
a periodic transversal displacement of the spicule axis at one 
particular height. 

If spicules are formed in thin magnetic flux tubes, then 
the periodic transverse displacement of the axis observed by 
Nikolsky & Platova <1971> probably was due to the propa- 
gation of kink waves. It is well known that transverse kink 
waves can be generated in photospheric magnetic tubes by 
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buffeting of granular motions (Roberts 119791 HoUweg I198U 
Spruit II 98 fl Hasan & Kalkofen 1999). Kink waves cause the 
displacement of the tube axis, therefore their propagation can 
be traced either by direct observation of the tube displacement 
along the limb as in Nikolsky & Platova (1971") or spectroscop- 
ically by the Doppler shift of spectral lines. The later possibility 
arises when the velocity of kink waves is polarized in the plane 
of observation. The periodic spatial distribution of Doppler ve- 
locities in spicules has been detected almost 20 years ago by 
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Fig. 1. The Doppler velocity spatial distributions for one of 
the height series is shown. Marked dots indicate the observed 
heights. The Doppler velocity has clear periodic spatial distri- 
bution, which indicates wave propagation. 

Khutsishvili J1986I I. Unfortunately, neither observers nor the- 
orists paid attention to the phenomenon which can be simply 
explained in terms of kink waves. In this letter, we reanalyze 
the height series of Ha spectra obtained by Khutsishvili ( 1986 ) 
and then model the observation as propagating kink waves. 

2. Observation 

The big (53 cm) coronagraph and universal spectrograph of 
Abastumani Astrophysical Observatory has been used to obtain 
chromospheric Ha line spectra at different heights (8 heights) 
above the photosphere (Khutsishvili 1986). Instrumental spec- 
tral resolution is 0.04 A and dispersion 1 A/mm in Ha. 
Observation has been caiTied out at the solar limb as height 
series beginning at 3800 km above the photosphere with a step 
size of 1 arcsec. Thus the distance ~ 3800-8700 km above the 
photosphere was covered. The exposure time was 0.4 s at lower 
heights and 0.8 s at higher ones. The duration of each height 
series was 7 s, with the total duration of the observation be- 
ing 44 min. More details about the observation can be found in 
KhutsishviU (.1986 ). 

We analyzed spatial distributions of Doppler velocities in 
selected He height series. Nearly 20% of measured height se- 
ries showed a periodic spatial distributions in the Doppler ve- 
locities. The Doppler velocity spatial distributions for one of 
the height series is shown in Fig.l. Periodic spatial distribu- 
tion is clearly seen, which indicates a wave propagation. The 
wave length can be estimated as ~ 3500 km. The maximal time 
difference between the observations at the lowest and high- 
est heights in one series is 7 s. If the wave propagates with 
the Alfven speed, being say ~ 50 km/s in the chromosphere, 
then during 7 s it will pass a distance ~ 350 km. So during 
one height series the wave will propagate along the distance 
which is less than the distance between consequent observed 
heights. Therefore the spatial structure of the wave velocity in 
one height series can be considered as nearly simultaneous. 

Thus the observations show the evidence of transversal 
wave propagation in spicules with typical wave length of ~ 
3500 km. The waves can be either kink or linearly polarised 




Fig. 2. Schematic picture of propagating kink waves in 
spicules. Due to the kink waves, the observed spectral line 
is Doppler shifted when the wave velocity is polarized in the 
plane of observation. 

Alfven waves (torsional Alfven waves in magnetic tubes will 
cause periodic broadening of spectral lines not the Doppler 
shift, see Zag arashvih 120031 . But as spicules are highly struc- 
tured phenomena (their width are ~ 1"), the observed Doppler 
shift can be caused due to the propagation of kink waves. 

Kink waves propagating along the magnetic tube lead to the 
oscillation of the tube's axis (see Fig. 2). Therefore if the veloc- 
ity of the kink wave is polarized in the plane of observation, 
then it results in the Doppler shift of the observed spectral line. 
The Doppler shift will have a periodic behaviour in height: at 
the same time, the spectral line will have blue shift at antinodes 
where the tube moves towards the observer and red shift where 
the tube moves in the opposite direction; at the velocity nodes 
the Doppler shift tends to zero. This means that kink waves 
manifest itself in the same spatial behaviour of Doppler shift 
as revealed by our height series (Fig. 1). However if the spicule 
axis is inclined with respect to the local vertical direction then 
a steady flow inside the spicule will shift the Doppler velocity 
at all heights with the same value. Indeed, some spicules show 
this behaviour. Three consecutive height series of Doppler ve- 
locity in one spicule are shown in Fig. 3. The velocity is shifted 
by ~ 12 km/s at all heights, although the wave signature is still 
seen. If the spicule is inclined by approximately 35 ° (Trujillo 
Bueno et al. '20051 then the real steady flow velocity will be ~ 
22 km/s, similar to the typical mass raising speed in spicules. 
Also Fig. 3 shows that the maximum of the Doppler velocity 
moves up in consecutive height series. This may indicate a 
wave phase propagation. The phase is displaced at ~ 1500 km 
in about 18 s giving the phase speed of ~ 80 km/s, very similar 
to the expected kink speed at these heights. 

3. The model 

The mechanism of spicule formation is still not well under- 
stood (see the recent review of Sterling 120001 and references 
therein). Therefore this letter does not address the spicule for- 
mation mechanism. Spicule life time is ~ 10-15 min, while the 
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Fig. 3. Three consecutive height series of Doppler velocity in 
one spicule. The time difference between the consecutive plots 
is ~ 8 s. The maximum of Doppler velocity moves up in con- 
secutive height series, which probably indicates a wave propa- 
gation. 

period of observed waves is much shorter: 70 s for the wave 
with phase speed of 50 km/s and wave length of 3500 km. In 
2-3 periods the wave will propagate along the whole spicule 
length. During this short time intervals spicules can be con- 
sidered as existing stable structures. Therefore the waves may 
propagate independently of the spicule formation mechanism. 

During short time intervals spicules can be modeled as thin 
magnetic flux tubes embedded in a field free environment, an- 
chored in the photosphere and persisted towards the corona. 
Propagation of kink waves in vertical magnetic tube embedded 
in a field free environment is governed by the equation (Roberts 
l2004l 



5^ 
df 



= Ct 



!^ ^ PO -Pe d£, 

po + Pedz 
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where ^ is the transverse displacement, po(z), pdz) are the 
plasma densities inside and outside the tube accordingly, g is 
the gravitational acceleration and c^. = Vpo/Po + Pe the kink 
speed. Here = Bo/ V^^Po is the Alfven speed with Bo(z) 
as the tube magnetic field. Note that this equation does not in- 
clude the flow along the tube, while spicules show continuous 
upward mass motion. Therefore equation (1) may correctly de- 
scribe the waves only in the co-moving frame (with the steady 
flow), and thus the wave phase speed will be Doppler shifted 
due to the mass motion inside the spicule, which may slightly 
alter the theoretical results. 

Spicule density and magnetic field show almost no spa- 
tial variation at observation heights (~ 3800-8700 km above 
the photosphere) and the stratification also can be neglected 
in this part of the solar atmosphere, then equation (1) gives 
the approximate dispersion relation c\k^^ - uP', where is the 
vertical wave number and o) is the frequency of kink waves. 
From the dispersion relation we may estimate the period of 
kink waves using the observed wave length and typical kink 
speed. The observed wave length of kink waves is of order 
~ 3500 km. Then for the kink speed of ~ 50 - 100 km/s at 
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Fig. 4. Comparison of observed Doppler velocity plot (Fig.l) 
and the numerical solution of the kink wave equation. Vertical 
axis shows the height in km and the horizontal axis shows the 
kink wave velocity in km/s. It is seen that the wave length of 
the computed kink waves at the photospheric level is ~ 1000 
km, which is comparable to the granular diameter. 



heights of 3800-8700 km above the photospheric, the period of 
kink waves can be estimated as ~ 35-70 s. 

It can be suggested that the source of observed kink waves 
resides in the lower part of the solar atmosphere. The cut-off 
period of kink waves due to stratification at the photosphere 
(for the pressure scale height of 125 km and plasma /3 of ~ 1) 
is ~ 660 s. So the expected period of kink waves is well be- 
low the cut-ofF value. Thus the kink waves with periods of ~ 
35-70 s may easily propagate upwards. As estimated frequency 
of the observed waves is much higher than the frequency of 
5-minute oscillations, then the oscillations can be ruled out to 
be the source of kink waves. Therefore, the only source which 
may excite the kink waves in the photosphere is the granula- 
tion. The photospheric granulation has been often suggested 
as the source of kink waves in thin magnetic tubes (Roberts 
\W7^ Spruit TWT HollwegTW Hasan & Kalkofen llg^ . 
Photospheric granulation is very dynamic even during the life 
time of one granular cell. If the magnetic tube, in which the 
spicule is formed, is anchored in the photosphere, then any per- 
turbation of the granular cell probably excites the kink waves 
with wave length similar to the cell diameter (see Fig.2). But 
the observed wave length (~ 3500 km) of the kink waves in the 
higher atmosphere is a few times longer than the mean granu- 
lar diameter, which is of order ~ 800 -1000 km. The discrep- 
ancy can be resolved by an increase of wave phase speed with 
height. Indeed, the kink speed being ~ 10 km/s at the photo- 
sphere increases up to ~ 50 km/s (or even more) in the higher 
atmosphere. Therefore it can be suggested that granular cells 
generate kink waves with a wave length comparable to their di- 
ameter, but the wave length increases due to increasing phase 
speed when waves propagate upwards. 

In order to show this, one should solve the kink wave equa- 
tion (1) along the whole solar atmosphere from the photo- 
sphere to the present observed heights. As equation (1) is ho- 
mogeneous with time, we may perform a Fourier expansion in 
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time in the form ^ = ^i(z) exp(/wf)i where w is the wave fre- 
quency. Then its numerical solution is straightforward; we only 
need to fix the spatial variations of the unperturbed parameters. 
Unfortunately, the height variation of the physical quantities 
is not well known. There are only suggested values of density 
and magnetic field in spicules and they are also controversial. 
Also the difference between plasma densities inside and out- 
side of spicule are not well determined. At the photospheric 
level the thermal equilibrium inside and outside the tube i.e. 
To(z) - Te(z) (here Tf,(z), Te{z) are plasma temperatures in- 
side and outside of tube) is a good approximation, which gives 
a higher density outside the tube than inside (in order to hold 
the transverse pressure balance at the tube boundaries). But in 
the higher atmosphere, say at heights > 2000 km, the plasma 
temperature is much lower in spicules than in the environment. 
On the-other-hand, the plasma density seems to be much higher 
in spicules than the surroundings. Due to this uncertainty in the 
medium parameters, it is better to give the height dependence 
of the Alfven speed and the density ratio inside and outside the 
spicule po/pe- Therefore we take the height dependence of the 
Alfven speed as 10 km/s at the photosphere and ~ 100 km/s at 
6000 km. For the density ratio we take 1/3 at the photosphere 
(being higher outside the tube), which increases up to 10^ at a 
height of 6000 km. With this spatial dependence of the Alfven 
speed and the density ratio we solve numerically the kink wave 
equation for different values of frequency a>. Then we fit the 
numerical solution to the observed curve of the Doppler ve- 
locity. The numerical solution is best fitted to the observation 
(Fig.4) when the wave period is ~ 40 s. From Fig. 4 it is seen 
that the wave length of the kink waves at the photospheric level 
is ~ 1000 km, which is comparable to the granular diameter. It 
clearly indicates a granular origin for the waves. 

4. Discussion 

Kink oscillations of coronal loops have been frequently ob- 
served in the solar corona by TRACE as periodic loop dis- 
placement in space. However in the lower atmosphere the ob- 
servation of kink waves is complicated due to their nearly in- 
compressible character. Here we report the spectroscopic ob- 
servation of kink waves in solar limb spicules. We have shown 
that nearly 20% of measured Ha height series at the solar limb 
(~ 3800 - 8700 km distance above the photosphere) show the 
periodic spatial distribution of Doppler velocities in spicules, 
which can be caused by propagating kink waves. The wave 
length at these heights is ~ 3500 km, which goes to ~ 1000 km 
at the photosphere if the height dependence of the kink speed 
is taken into account. The estimated wave period is ~ 35-70 s. 
The computed wave length at the photospheric level is com- 
parable to the granular size, therefore the granulation is proba- 
bly the most plausible source for the wave excitation. Then the 
granulation will excite kink waves in any thin magnetic flux 
tube anchored in the photosphere. As we already noted only 
~ 20% of measured height series show periodic spatial distri- 
bution of Doppler velocities. Indeed spectroscopic detection of 
kink waves is possible only when the wave velocity is polarised 
in the plane of observation. But the granulation will excite the 
waves with any direction of polarisation; there is no preferred 



direction. Therefore 20% of measured height series is a good 
percentage for the observation of waves. 

It must be noted that some spicules show a multi- 
component structure (i.e. several spicules are located close to- 
gether) mostly at lower heights, which may give the impression 
of a velocity shift (Xia et al. 2005). Therefore we choose only 
the spicules with well defined single-component structure. We 
also note that the data analysis has been made only for selected 
height series. Now we are doing the data reduction of the whole 
observation, which will be analysed in the near future. 

Here we do not address the mechanism of wave excitation 
in detail as it is beyond the scope of this letter. But we sug- 
gest that granular cells will frequently excite the kink waves 
with the wave length comparable to their diameter in anchored 
magnetic tubes. The waves may propagate upwards carrying 
the photospheric energy into the corona, thus may cause sig- 
nificant input into the coronal heating. Therefore detailed spec- 
troscopic search of waves with period of ~ 40-60 s and wave 
length of ~ 3000-4000 km is desirable to be performed in the 
future. 
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